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Abstract 



In this paper we derive a general expression for the acoustic Casimir pressure between two 
parallel slabs made of arbitrary materials and whose acoustic reflection coefficients are not 
equal. The formalism is based on the calculation of the local density of modes using a 
Green's function approach. The results for the Casimir acoustic pressure are generalized to 
a sphere/plate configuration using the proximity theorem, 
pacs: 43.25. Qp, 43.10.Ln 
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I. INTRODUCTION 



Due to quantum vacuum fluctuations, two para^ 
This phenomenon is known as the Casimir force 



lei neutral plates will attract each other. 
Although a small force, it has been 
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measured accuratelyusing torsion balances, atomic force microscopes and micro mechanical 



oscillators 



BflQ 



|. An acoustic analog to the Casimir effect was reported in 1996 by 
Larraza and collaborators 13? 5 where two parallel plates, placed at a distance L, were 
subjected to a broadband noise back ground. The plates were observed both to attract and 
to repel each other, depending on the separation distance and the broadband noise cutoff 
frequencies. Following Casimir's method 1], a theory for the acoustic Casimir force was 
developed assuming perfectly reflective plates with approximations that turned out to be 
valid for the frequency range, material and plate thickness that were used in the experiment. 
In this work we derive a general expression for the Acoustic Casimir pressure for materials 
with arbitrary impedances by calculating the density of modes between the plates using 
the Green's function formalism borrowed from the electromagnetic case. We also present 
an acoustic analog to the proximity theorem to calculate the Casimir pressure between a 
sphere and a plate. 

II. THEORY 

Consider two different parallel slabs labeled i — 1, 2 of thickness d\ t 2, separated by a 
distance L along the z axis. The slabs are parallel to the x — y plane and have an arbitrary 
acoustic reflectivity ri j2 (Fig. la). 

For a perfect acoustic reflector, the radiation pressure of a wave of intensity I and speed 
c impinging on the slab is given by (J 0| 

P=^cos 2 (#), (1) 

where 9 is the angle of incidence. 

As in Ref. jf], Q, 0] we consider broadband acoustic noise of constant spectral intensity 
I u in the frequency interval [ui, u 2 ], and its spectral representation in the wave vector space 
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where k 2 = (uj/c) 2 = k 2 + k 2 + k 2 . 

The total radiation pressure on a plate for non perfect reflectors, using Eq. (JTJ) is 

I* f k 2 
Po = / dk x dkydk z -j-£. (3) 

Between the plates, the total radiation pressureP in is determined by the allowed modes 
that satisfy the boundary conditions at the plate surfaces. If the plates are large enough, 
k x and k y take on continuous values. For perfect reflectors the normal component of the 
wave vector takes the values k z = nir/L, where n is an integer. Thus, the calculation of the 
energy density is reduced to integration over kx and ky and summation over n. The Casimir 
pressure results from the difference P{ n — P a . 

For arbitrary materials the mode summation is no longer direct, since we are no longer 
allowed to specify Dirichlet boundary conditions that restrict the allowed modes, and it 
becomes necessary to calculate the total density of modes T>(k z ) and integrate over all wave 
vector space. To do this we use a Green's function approach. 

The wave equation for the velocity potential can be written as the eigenvalue equation 
—d 2 (j) = k 2 (p, with eigenvalues k 2 . Let A n be an eigenvalue for the eigenfunction <f> n . In 
terms of the velocity potential 0, the particle velocity v z and the fluid density p and with 
the definition s v r = d z <f) n and p = —pdt<f) n we can write the normal component of the wave 

y. 

w n = P - ((d z <p n ) 2 + k 2 z <&) . (4) 



stress tensor 10 as 



The total contribution to the stress tensor is obtained by adding up all modes and inte- 
grating over all possible values of k 2 . 



P 

W= 2 



J dk 2 z S(k 2 z - A n ) {{d z <p n f + k 2 z <p 2 n ) (5) 



where we have assumed an harmonic behavior of the potential and the Dirac's delta 
function is introduced since only the eigenmodes contribute to w. Now, using the identity 

' Ptt -V -^(k 2 z -X n ), (6) 



k+ 2 - \ n k\ - \ n 



with k^ 2 = \im v ^ Q (k 2 + ir/) we can write Eq.© as 



4 



In this equation, we can identify the spectral representation (or eigenfunction expansion) of 



the Green's function and its derivative 



and we interpret the quantity 



V k 2 = --Im(G(z, z) + d z d z G(z, z)), (8) 

7T 

as the density of modes. Another way of understanding the result is as follows. The basic 
definition of denisty of modes in terms of the Green's funtion is obtained from Eq. (JHJ) as 
—ImG{z ) z)/7i. The acoustic pressure obeys the wave equation and with th eappropriate 
boundary conditions, we can obtain the Green's function G p and thus the density of modes. 
From the acoustic stress tensor component (Eq.(@J)), besides the pressure field there is a 
contribution from the velocity field. Let this field have an associated Green's function G v 
. The total density of modes of the system —Im(G p + G v )/n of both fields. Writing the 
pressure and velocity in terms of the scalar potential yields Eq. (JHJ). This is equivalent to 
what happens in zero point Casimir effect where the density of modes comes from adding 
the contribution of the electric field plus that due to the magnetic fields, and both fields are 
related through a constitutive equation (Maxwell's equations). 

To construct the Green's function for the velocity potential we can use the standard 
definition 

G k 2(z,z) = — , (9) 

where W is the Wronskian and 

(t)<(z) = e~ ikzZ + r,e ,fe2 

(10) 



(j)>( z ) = e ifc ^- L ) + r 2 e- ikz{z ~ L \ 

are the solutions to the one dimenional wave equation where the super-index (<, >)represents 
the smaller and larger of z and z' respectively. 

Substitution of the potentials (Eq. (JTOJ)) into Eq.(JHJ) and Eq. (JHJ) yields the local density 
of modes 



where we have obviated the dependence of the reflectivities with wave vector. The density 
of states Eq. fjllJ) was obtained from the Green's function for the Helmholtz equation with 
eigenvalues k 2 z . However, we are interested in the density of states for k z . This is simply 
obtained from T>^2 = d(k z )T>k z = 2k z T>^, or 

1 + nr 2 e 2ikzL ] 

(12) 



V kz = -Re 

7T 



The radiation pressure due to the inside modes can now be written as 

Pin = ^r / dkxdk y dkzT>k ^- ( 13 ) 

In the limit of perfect reflectors r — > 1 the density of states becomes ^-5{k z — nko) where 
k Q = n/L. Thus, from Eq. (JT3j) the pressure due to all the modes is 

k I f 

Pin = ~^T / dk x dk y dk z ^ 5(k z - nk )-^, (14) 



or 

Pm = ~~2^r S dkxdky ^ W+W+^W' (15) 

that is the same as that obtained by Larraza and collaborators ((j, 3>Q|. Finally, the acoustic 
Casimir force per unit area / = Pj n — P out takes the form 

f = (/ dk - dk v dk ^JZ^) > ( 16 ) 

where ^ = (r!r 2 exp (2ik z L))~ l . Notice that it is enough to know the separation between 
the slabs and the reflectivities to determine the acoustic Casimir force. 

To illustrate the application of Eq. (jlfjjl . in Fig. 2 we plot the force vs separation for 
two identical slabs with constant reflectivities r = 1, .8, .7. In all cases the force goes from 
attractive to repulsive as the separation increases. The magnitude of the force is not only 
related to the reflectivity but also to the finite bandwidth being used. If the band width 
extends from zero to infinity, the acoustic Casimir pressure for a perfect reflector —ttI^/AL is 
always attractive. From Eq. (fT7)|) if we integrate over all frequencies the force is also always 
attractive and as the reflectivity decreases the force does too 3] in all cases. Without loss of 
generality we have assumed a constant value of r within the bandwidth under consideration. 
However, the formalism is valid even when the reflectivity shows a strong dependence with 
: requency. The bandwidth and intensity used in these calculations are the same as in Laraza 
f|. Interestingly, even if we consider a finite frquency bandwidth it is possible to obtain a 
purely attractive force if we consider the force between a surface with reflectivity r\ = 1 and 
a pressure release surface r 2 = — 1. In this limit the force is constant and always attractive 
since Dj tz 0, as can be seen from Eq. ()12|) and the external pressure field pushes the 
plates together. 
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III. THE PROXIMITY THEOREM IN ACOUSTICS 



Practical measurements of the Casimir force due to zero point energy fluctuations are 
done between a large sphere and a plane due to the difficulty of keeping two plates parallel 
at the submicron scale . The force between a large sphere and a plane (see Fig. lb) 

is calculated using the proximity theorem[16] or Derjaguin approximation 17], which states 
that 

the force between two smooth surfaces as a function of the separation degree of 
freedom is proportional to the interaction potential per unit area £ between 
two flat surfaces, the proportionality factor being 2tt times the reciprocal of 
the square root of the Gasussian cuadrature of the gap width function at the 
point of closest approach . 

For a sphere-plane system, the force F sp is obtained from the Casimir free energy per unit 
area between two parallel plates £ as 

F sp = 2tcR£, (17) 

where R is the radius of the sphere. The proximity theorem is valid provided L/R < 1, being 

lough the limit of L — > can be described 



121 ] . is that there are no bounds on how big 



L the closest distance between the surface, alt 
by the proximity theorem. A current problem 
L/R has to be in order to obtain the correct result. Experimentally this becomes difficult 
at the submicron scale. The acoustic analog of the Casimir force provides a simpler (not 
necessarily easier) way of solving this problem, since as shown by Laraza [6j the scale of the 
acoustic experiments allows a more precise control of the involved parameters. Furthermore, 
the proximity theorem is valid for any interaction. In the acoustic case, the free energy per 
unit area for parallel plates is 

S = {^J dKdkydk^Re (Zn(l - r x r 2 e 2 ^ L ))) . (18) 

This expression for £ is such that the force (Eq. ()16|0 is given by / = —d£/dL. Thus, the 
force between a sphere and a plane is 

fp° = ^rj dk x dk y dk z ^Im (ln(l - r ir2 e 2 ^ L ))) . (19) 
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To show the application of this approximation, in Fig. 3 we have plotted the force between 
two parallel plates and the force between a 20cm sphere and a plate. We observe that 
although the proximity theorem gives a correct behavior and overall oreder of magnitude 
for the force, the region in the limit of L approaching zero is not well described. 



IV. CONCLUSIONS 



We have derived a general expression for the Acoustic Casimir force between two parallel 
slabs with arbitrary acoustic properties characterized by the reflection coefficients of the ma- 
terial. We also extended our results to include the force between a sphere and a plane. The 
expression we obtain for the Casimir force is convenient for calculations since it depends 
mainly on the reflection coefficients that can be obtained straightforwardly in theoretical 
computations or through experimental work. Our approach is analogous to the electromag- 
netic dielectric case, so this formalism is equivalent to the Liftshitz formula In the 
limit of a perfectly reflective plate our results agree with those of L . This formal- 
ism can be extended to the case of highly porous materials or viscous propagation media, 
although the calculations involved can be of increasing difficulty. It must be pointed out 
that for the case when the material is deformed by the wave, this density of states approach 
is no longer valid: the reflection coefficient is angular dependent, and since the angle itself 
is time dependent the use of a static density of states would be incorrect. Also, we have 
excluded the possible effects of roughness. 

The crucial difference between the system we consider and the original treatment [f| is 
the inclusion of the density of states through the Green's function method. The analytical 
interpretation of the density function gives a deeper insight into what really happens in a 
non perfect reflector. For a perfect reflector the density of modes consists of a series of 
Dirac's deltas. As the reflectivity decreases from unity, the resonance bands increase in 
width which is heuristically equivalent to a spatial diffusion of the nodes that appear inside 
the resonant cavities. 

Although the use of perfectly reflectiving plates is a good approximation in some ex- 
perimental situations, this is not the case for other bodies (such as rubber, as an extreme 
example) hence our efforts to broaden the horizon of application. As an example, we have 
consider the possibility of using acoustic experiments to prove the validity of the proximity 
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theorem. Additionally, this treatment could allow for a larger range of experimental vs. 
theoretical comparison in this field where, as noted by Larraza et al, the possibility of direct 
technological application of the acoustic Casimir effect is considerable. 
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FIG. 1: a) Geometry and coordinate system for the two parallel plate configuration and b) for the 
sphere-plane configuration. 

FIG. 2: Acoustic Casimir force between two parallel plates for different values of the reflectivity 
r assuming both plates are equal. The values of r are indicated in the figure. The intensity and 
bandwidth are the same as in the experiments of Laraza [(]]). 

FIG. 3: Acoustic Casimir force between a sphere (R=0.2 m) and a plane for reflectivity r = 1 
(dashed line). For comparison the force between two parallel plates is also shown (solid line). 
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